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We experimentally study the effect of near field coupling on the transmission of light in ter¬ 
ahertz metasurfaces, possessing slightly distinctive SRR resonances. Our results show that the 
interplay between the strengths of electric and magnetic dipoles, modulates the amplitude of result¬ 
ing electromagneticahy induced transmission, probed under different types of asymmetries in the 
coupled system. We employ a two-particle model to theoretically study the inhuence of the near 
field coupling between bright and quasi-dark modes on the transmission properties of the coupled 
system and we hnd an excellent agreement with our observed results. Adding to the enhanced 
transmission characteristics, our results provide a deeper insight into the metamaterial analogues 
of atomic electromagneticahy induced transparency and offer an approach to engineer slow light 
devices, broadband hlters and attenuators at terahertz frequencies. 


Light-matter interaction has been a subject of intense 
research over past several decades, since it allows to probe 
the resonance and the off-resonance properties of the ma¬ 
terials over large part of the electromagnetic spectrum. 
Until late twentieth century, light-matter interaction in 
the terahertz part of the electromagnetic spectrum was 
the least explored. With the advent of metamaterials [1]- 
[3], which exhibit structure dependent resonance proper¬ 
ties, have become excellent candidates for probing such 
resonant and off-resonant interactions at terahertz fre¬ 
quencies. Metamaterials are composed of periodic ar¬ 
ray of sub wavelength sized met a-atoms, which exhibit 
strong near-field coupling that can carry the interaction 
energy over to the far field regimes. Super lens [H |5], 
hybridization[6H9], Fano-couphng[T0HI2] and the classi¬ 
cal analogue of electromagnetically induced transparency 
fEITj pRUTT] have been studied and demonstrated using 
the near field coupling within the metamaterials. Re¬ 
cently, there have been a enormous interest in the near¬ 
field coupling in terahertz metamaterials, which show 
EIT like transmission p!8l - [2Q] and ultra high Q Eano 
resonances [ 2 TI [22] . which find significant applications in 
the terahertz sensing [23[ |24] and broadband communica¬ 
tion technologies [25]. 

Electromagnetically induced transparency is a quan¬ 
tum interference effect, which was first observed [26] in a 
three level atomic system, owing to the destructive inter¬ 
ference between the possible excitation pathways. Later 
its analogue was extended to the classical systems [27] . 
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since then EIT effects have been observed in various clas¬ 
sical systems, including metamaterials. photonic 
crystals [28]. micro ring resonators [29l [30] and all dielec¬ 
tric metasurfaces [3T] . There have been a few reports on 
tailoring the classical analogue of EIT using metamate¬ 
rials at microwave [32ll34] . terahertz [35II37] and optical 
frequencies [T3I [38] . either by tuning the near field cou¬ 
pling or by changing the material properties. Manipula¬ 
tion of EIT in classical systems will allow us to precisely 
tailor the group velocity [T9l|38] and the delay bandwidth 
product [32] of the transmitted pulse. Moreover, it pro¬ 
vides a clear picture of the coupling mechanisms in the 
classical analogue of EIT, that can help us in drawing 
the closest analogy between the classical and the quan¬ 
tum systems. 

In this letter, we experimentally demonstrate the en¬ 
hancement and suppression in transmission of the fields 
at terahertz frequencies, by manipulating the near field 
coupling between the radiative dipole ring and the sub¬ 
radiant quasi-dark split ring resonator (SRR) ring in 
metasurfaces under different type of asymmetries of 
the metamolecule. Systems with enhanced transmis¬ 
sion show a considerable increase in the delay bandwidth 
product (DBF) at the transmission peak. We quanti¬ 
tatively interpret underlying phenomenon by using the 
two-particle model, which shows an excellent agreement 
with the observed experimental results and is discussed 
in detail in the following sections. 

Metamaterial unit cells (shown in Eig.l), consist of a 
metallic split ring resonator (SRR) surrounded by a con¬ 
centric metallic closed square ring resonator (CRR), both 
having a thickness of 200nm. Samples were fabricated us¬ 
ing photolithography technique, where 200nm thin layer 
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FIG. 1. (color online). Metasurfaces showing (a) symmetry (MS- 
1), (b) SRR-position asymmetry (MS-2), (c) SRR-gap asymmetry 
(MS-3) and (d) SRR-position and SRR-gap asymmetry (MS-4), 
are graphically displayed. In the case of (b) and (d) SRR ring is 
displaced upwards by 4/im from the symmetric axis (x), whereas in 
(c) and (d) SRR gap is displaced sidewards by 5p,m from the sym¬ 
metric axis (Y). All four metasurface samples have same material 
dimension: L, dO/xm; S, 20/rm; g, 4/rm; w, S/xm; Periodicity (Px) 
of the unit cell is 50p,m with substrate thickness {Pz) of 640/i,m. 


of aluminium is deposited on 640/im thick n-type silicon 
substrate (e = 11.68). Structural symmetry of the meta¬ 
molecule unit cells were broken to study the impact of the 
near field coupling in the asymmetric metasurface array. 
Metasurfaces, MS-3 and MS-4, are the broken SRR-gap 
symmetry structures of MS-1 and MS-2 respectively (cat¬ 
egorised as SRR-gap asymmetry), whereas, MS-2 and 
MS-4 are the broken SRR-position symmetries of MS- 

1 and MS-3 respectively (categorized as SRR-position 
asymmetry). For the metasurfaces (MS-3 and MS-4) 
with SRR-gap asymmetry, the capacitive split gap in the 
SRRs of MS-1 and MS-2 is displaced horizontally (along 
the x-axis) by 5/rm from the ^/-symmetry axis, whereas 
for the metasurfaces with SRR-position asymmetry (MS- 

2 and MS-4), position of the inner SRR ring in MS-1 and 
MS-3 is displaced upwards (along the ^-axis) by 4/im 
from the x-symmetry axis (see Fig.l). The design of 
each metamolecule was chosen such that the fundamen¬ 
tal resonance frequencies of their constituent resonators, 
exhibiting highly contrasting resonance linewidths, fall at 
the same frequency, which is essential to realize the EIT 
like behavior in classical systems. Fig.2(a) depicts the 
contrasting resonance linewidths for the CRR ring and 
the SRRs. The measured Q-factor of CRR is 1.2 which is 
an order of magnitude lower than that of the Q-factor of 
inner SRR. For the symmetric SRR, the Q-factor is 10.6, 
whereas the asymmetric SRR (SRR with displaced gap) 
has a Q-factor of 11.3. Here, we would like to stress that 
both the resonators, CRR and the SRR interact with the 
incoming electric field {E^) but with very different cou¬ 



FIG. 2. (color online), (a) Transmission vs frequency graph dis¬ 
playing the individual resonance dips for the outer closed resonator 
ring (CRR) (black), the inner symmetric SRR (red) and the inner 
asymmetric SRR (blue) rings. Electric field strengths for the sym¬ 
metric (b) and the asymmetric (c) SRR rings are also shown. 


pling strengths, which is determined by the Q-factors of 
their fundamental resonances. The SRR ring with sharp 
LC resonance (higher quality factor) is termed as ‘quasi¬ 
dark mode’, because of its weak coupling to the incoming 
light field (E^,), where as the CRR ring with broad dipo¬ 
lar resonance (lower quality factor) behaves like a ‘bright 
mode’ that couples strongly to the incoming light field 
{E,). 

Figures 3(a)-(d) show the sharp transmission spec¬ 
tra for each of the metamaterial samples. Spectra 
were recorded using 8f confocal terahertz time domain 
spectroscopy (THz-TDS) for the incident electric field 
{Ex)^ polarized along the gap (x-axis) of the SRR ring. 
Recorded transmission time domain signals were con¬ 
verted to frequency domain data using FFT and nor¬ 
malized to the transmission of the bare silicon substrate 
(as reference). Such narrow transmission peaks result 
from the destructive interference between the classical 
resonators (CRR and SRR) due to their strong near field 
coupling. Corresponding numerical simulations were 
carried out using the commercially available CST MI¬ 
CROWAVE STUDIO Maxwell equation solver and the 
data matched well with our measured results (insets of 
Fig.3). 

Eig.3(a) depicts EIT like transmission spectrum for 
MS-1, where, between two transmission dips (at 0.83THz 
and 1.02THz) a sharp transmission peak is observed at 
the frequency 0.088THz, which signifies the hybridization 
model of the plasmonic coupling via near fields of the in¬ 
dividual resonators. The resonance dip at the lower fre¬ 
quency is the (i) antisymmetric mode and the higher fre¬ 
quency resonance dip is the (iii) symmetric mode. This 
behaviour is verified by simulating the surface current 
distribution as shown in the Eig.4(a) for MS-1, where 
at the transmission dips ((i) & (hi)), surface currents 
in the SRR and CRR run antiparallel and parallel to 
each other respectively for the (i) antisymmetric and the 
(iii) symmetric modes. Existence of the antisymmetric 
mode at the lower frequency signifies strong transverse 
dipole-dipole interaction within the coupled system. This 
dipole-dipole interaction can be electric as well as mag- 
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FIG. 3. (color online). Terahertz transmission spectra for the 
incident Ex field. Figures (a), (b), (c) and (d) show the experi¬ 
mentally measured as well as simulated transmission (inset) curves 
for metasurface samples MS-1, MS-2, MS-3 and MS-4 respectively. 
Experimentally measured transmission curves are fitted with corre¬ 
sponding the analytic model data (colored magenta), (i), (ii) and 
(iii) in (a), simultaneously represents the antisymmetric, transmis¬ 
sion and symmetric modes of the plasmonic hybridization, which 
is detailed in Fig.4(a). 


netic as the respective dipoles (electric/magnetic) are 
aligned in the transverse way to the mode separation 
line. Electric dipoles are aligned antiparellel to each other 
in the plane of the material and the induced magnetic 
dipoles are directed perpendicular to the metamaterial 
surface. 

The observed transparency peaks (Fig.3) can be well 
explained using the induced surface currents and elec¬ 
tric field distribution as shown in Fig.4(a & b). Fig.4(a) 
shows the the induced surface currents in the coupled 
system, where at the (ii) transparency peak induced cir¬ 
cular currents in the SRR destructively interferes with 
the dipolar currents of CRR. As an evidence of the de¬ 
structive interference, we observe extremely weak surface 
currents in the CRR that are opposite in direction to 
the surface current in the SRR. The same can be ex¬ 
plained using electric field distributions (see Fig.4(b)) 
at the (i)antisymmetric resonance, (ii)transparency peak, 
and the (iii) symmetric resonance mode for each metasur¬ 
face sample. By comparing the maximum electric field 
strength in the SRR gap at the (ii) transmission peaks 
with the transmission dips in (i) and (iii), clearly suggests 
that within the transparency peak SRR’s LC resonance 
dominates the dipolar CRR resonance. Thus the sur¬ 
face current and the electric field appears to be localized 
within the SRR of the coupled system. 

Introducing SRR-gap asymmetry in the system (MS-3) 
lead to a suppression in the transmission amplitude of the 
resulting transparency peak (Fig.3(c)). This suppression 
in the transparency peak is solely due to the weakened 


SRR resonance (ref blue curve in Fig.2(a)), because of its 
structural asymmetry. Displacing the gap to one end of 
the SRR arm results in rather a flaccid and weak electric 
field distribution (Fig.2(c)) at the SRR gap. This weak¬ 
ens the capacitive coupling within the SRR ring. As a re¬ 
sult, the effective strength of the quasi-dark mode (SRR) 
decreases and results in a reduced transmission. Owing 
to the interference effects in EIT phenomenon, this result 
can be seen as analogous to the waves’ interference, where 
decreasing the amplitude of one of the wave, results in 
the decreased strength of interference pattern. The same 
explanation holds true for the observed suppression in 
the transmission for MS-4 (Fig.3(d)) compared with the 
transmission of MS-2 (Fig.3(b)). 

Upon introducing SRR-position asymmetry in the sys¬ 
tem, transmission through MS-2(4) shows an enhance¬ 
ment over the transmission observed for MS-1(3), as 
shown in Fig.3(b & d). Enhanced transmission is due 
to increased coupling between the CRR and the inner 
SRR, which results in strong destructive interference of 
the two fields at the transmission peak. This effect can be 
explained using Fig.4(c), which depicts the change in the 
strength of induced surface currents at the transparency 
peaks, for all the four coupled metasuface structures. For 
MS-1(3), strongly confined fields of the SRR induces op¬ 
posing currents in CRR, which results in a destructive 
interference of the fields, giving rise to a sharp trans¬ 
mission peak. As the SRR ring is moved upwards (MS- 
2(4)), the separation between the bright (CRR) and the 
quasi-dark (SRR) modes decreases, which result in re¬ 
duced current density in CRR. This is because, when the 
SRR gap (stronger E-field confinement) comes closer to 
the bright mode (top arm of the CRR), electric coupling 
between the two modes dominates the interaction. This 
reduced separation enhances the effective coupling in the 
system that leads to a strong cancellation of the oppos¬ 
ing currents (enhanced destructive interference) within 
the coupled modes, resulting in an increased transmis¬ 
sion of the incident field. This effect is also reflected in 
the F-field distribution diagram as shown in Fig.4(b)(ii), 
where at the transparency peaks for the metasurfaces 
MS-2 and MS-4, the in the coupled system is decreased 
(compared with MS-1 and MS-2) as the result of the en¬ 
hanced destructive interference between the two modes. 
Thus the entire system behaves as superadiative system 
(smaller Q-value) due to the increased coupling within 
the coupled system. On the other hand, vice versa holds 
true, when the SRR ring is displaced downwards. The 
observed frequency red shift of the transmission peak and 
the asymmetric mode for the MS-2 and MS-4 (see Fig.3), 
indicates the increased electric field strength within the 
coupled metamaterial system [40]. Thus by moving the 
quasi-dark mode relative to the bright mode, we can 
modulate the electromagnetically induced transparency 
by tailoring the strength of the electric and magnetic cou¬ 
pling within the system. 

Effect of coupling on the transmission of light in the hy¬ 
bridized metamaterial system is theoretically studied by 
using two particle model. In the following coupled differ- 
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FIG. 4. (color online). (a) Induced sur¬ 
face currents in SRR and the CRR at 
transmission dips ((i), 0.83THz and (iii), 
1.02THz) and at the transmission peak 
((ii), 0.88THz) for the metasurface MS- 
1, are shown. 

(b) Fa;-field distributions at two trans¬ 
mission dips ((i) & (iii)) and at the trans¬ 
mission peak (ii) are shown for each meta¬ 
surface samples. 

(c) Depicts the surface current distribu¬ 
tions for the metasurfaces MS-1, MS-2, 
MS-3 and MS-4 at their respective trans¬ 
mission peaks. 


ential equations, we consider both particles (bright (x^) 
and quasi-dark (xd)) interact with the incoming electric 
field E = 


Xb{t) + + ‘xlxbit) + ^'^Xd{t) = — (1) 


Xd{t) + ldXd{t) + uJdXd{t) + n'^Xbit) = ( 2 ) 

Here, (Q, qd), {M, nid), {uib, uid) and ( 7 b, 7 ^) are the 
effective charge, effective mass, resonance angular fre¬ 
quencies and the loss factors of the bright and the quasi¬ 
dark respectively, ft defines the coupling strength be¬ 
tween the bright and quasi-dark particles. In the above 
coupled equations, we substitute Qd = ^ and rud = 
which gives the relative coupling of incoming radiation 
with the bright and the quasi-dark modes. Here, A and 
B are numerical constants. Now by expressing the dis¬ 
placements vectors for bright and quasi-dark modes as 
xi) = and Xd = CdC^^^^ we solve the above coupled 

equations ( 1 ) and ( 2 ) for and Xd- 


^ + iujjd)){QE/M) 

and 

^ + {B/A){uj^ - ul + iuj^b)){QE/M) 

The linear susceptibility (y), which relates the polar¬ 
ization (P) of the particle to the strength of incoming 
electric field (E) is expressed in terms of the displace¬ 


ment vectors as, y = ^ 

X = Xr + ^Xi 

K A{B^l)n^^A^{u:^-ul)^B{u;^-ul) 

A^B - ^6 + iujjb){^‘^ “ 

. A^Xd + Bjb X 

- ^6 + iujjb){^^ “ 

Here, Re[y] represents the dispersion and Im[y] gives 
the absorption (loss) within the medium. We fit 1- 
Im[y], to the experimental data shown in Fig.3 (colored 
magenta), which represents the transmission through a 
medium. For the fit, the values of the loss factors and 
jd are obtained from the linewidths of the curves shown 
in Fig.2(a), which are calculated to be around 3 x 10^^ 
rad/sec and 5 x 10^^ rad/sec respectively. The coupling 
strength ft for each transmission curves is calculated us¬ 
ing the formula given in Ref. [33], which can also be 
derived from Eqn.5 at the stop-band frequencies uj±^ for 
a loss less medium (assuming Ub^ood = Using Eqn.5, 
at the transparency peak {uj = ujt) where Re[y] = 0, we 
get , ujd ~ ujt for larger A. At the stop-bands where 
Im[y]= oc, by using Eqn.5 we can arrive at the expres¬ 
sion for uJb (assuming ood = cjt), (^b± = ^ ^ . By 

substituting the calculated values for 7 ^, 7 ^^, uJd, ft 
and by putting B = 2 (mass of SRR is half the mass of 
CRR) and iX = 4 x 10 ^^ (amplitude offset) in Eqn.5, we 
find an excellent fit for the transmission curves shown in 
Fig.3, for parameter A = 40. This suggests that, within 
these coupled metamaterial systems, interaction of the 
quasi-dark mode to the incoming radiation is nearly 20 
times smaller than that of the bright mode. This seems 
reasonable considering an order of magnitude difference 
in their respective Q-factors. Thus, susceptibility expres¬ 
sion given in Eqn.5 provides close to experimental situ- 
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Coupling strength, Cl (xio^i2 rad/sec) Frequency (THz) 


FIG. 5. (color online), (a) fit to the variation of Q-factors 

of the transmission curves with the coupling strength CL^ under var¬ 
ious SRR-position asymmetries of the metasurface samples. Brown 
curve represents Q-factor variation for metasurfaces with symmet¬ 
ric SRR ring, and green curve represents metasurfaces with asym¬ 
metric SRR ring, (b) Both experimental and simulated group delay 
data for the incident radiation within the transparency peak. 


ations of the present system under consideration with 
very small deviations. To further evaluate coupling ef¬ 
fects in these systems, we study the influence of cou¬ 
pling strength(f]) on the Q-factors of the transmission 
curves. We find that the Q-factors obtained from the 
simulated transmission curves for the samples with dif¬ 
ferent SRR positions (SRR-position asymmetries) follow 
^ variation (Fig.5(a)), as predicted by the two particle 
model [33]. 

Figure 5(b) shows a variation of the measured group 
delay values for the transmission curves given in Fig.3. 
Experimentally measured values for group delay, delay 
band-width product (DBF) and the respective Q-factors 
for the transmission curves for all the four metasurface 
samples are listed in Table. 1. From the data we see that 
as coupling strength increases, DBF increases and the 
corresponding Q-factor decreases and vice versa. For ex¬ 
ample, metasurface-2 that displays stronger mode cou¬ 
pling compared to all other samples, possesses maximum 
DBF and minimum Q-factor. In addition to the tun- 
ability of the group delay, the proposed planar slow light 


metasurfaces offer improved DBF at the terahertz fre¬ 
quencies, indicating enhanced buffering capabilities for 
the broadband telecommunication networks. 


Samples 

Group delay 

tg = ^(Ps) 

DBP 

{tg X A/) 

Q-factor 

Sample-1 

2.77 

0.21 

11.74 

Sample-2 

2.52 

0.217 

9.83 

Sample-3 

2.63 

0.174 

13.35 

Sample-4 

2.35 

0.18 
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TABLE I. Calculated group delay, DBP and Q-factor values for 
the experimentally measured transmission curves shown in Fig.3. 


In summary, we have experimentally studied the near 
field effects on the transmission properties of planar 
metasurfaces under two types asymmetries in the system. 
We demonstrated that, introducing SRR-gap asymmetry 
suppresses the transmission, whereas in the system with 
SRR-position asymmetry, where the SRR gap is moved 
towards the bright mode, enhanced electric coupling will 
enhance the transmission. Our results show that the 
transmission amplitude can be modulated by varying the 
strength of electric and the magnetic dipoles. The ef¬ 
fect of interplay between the strengths of electric and 
magnetic dipoles on the transmission provides a deeper 
insight into the classical analogy of the quantum interfer¬ 
ence effect. Using two particle model, we quantitatively 
analyzed the interactions within the system and the an¬ 
alytic data showed an excellent agreement with the ob¬ 
served results. The model also provides the relative cou¬ 
pling of the quasi-dark (SRR) mode to the incoming field 
in the near field coupled system. The proposed asym¬ 
metric planar slow light metasurfaces with tunable trans¬ 
parency characteristics will allow us to precisely control 
the group velocity of the pulse within the medium. They 
can be readily applied in broadband terahertz technolo¬ 
gies and show potential applications as variable power 
attenuators, broadband filters, and compact delay lines 
for the terahertz waves. 
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